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Abstract 

In this paper we study the problem of the nurnorical calculation (by Monte Carlo Methods) of the 
effective diffusivity for a particle moving in a periodic divergent-free velocity filed, in the limit 
of vanishing molecular diflfusion. In this limit traditional numerical methods typically fail, since 
they do not represent accurately the geometry of the underlying deterministic dynamics. We 
propose a stochastic splitting method that takes into account the volume preserving property 
of the equations motion in the absence of noise, and when inertial effects can be neglected. 
An extension of the method is then proposed for the cases where the noise has a non trivial 
time-correlation structure and when inertial effects cannot be neglected. Modified equations are 
used to perform backward error analysis. The new stochastic geometric integrators are shown 
to outperform standard Euler-based integrators. Various asymptotic limits of physical interest 
are investigated by means of numerical experiments, using the new integrators. 
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1. Introduction 

Understanding the transport properties of particles moving in fluid flows and sub- 
ject to molecular diffusion is a problem of great theoretical and practical importance 
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[7ll5j with applications in, for example, atmosphere, ocean science and chemical engi- 
neering f 5l31l21j . In the case where inertial effects can be neglected, the equation of 
motion for the particle is 

x = v{x,t)+aW. (1.1) 
Here a; € K'^ and is a standard d-dimensional Brownian motion, v{x,t) is the fluid 
velocity field (which we take to be incompressible) and a the molecular diffusivity. We 
will refer to this equation as the passive tracers model. 

It can be shown using multiscale/homogenization techniques |2l23j that when the ve- 
locity field v{x,t) is either periodic or random with sufficiently good mixing properties. 



the long time, large scale dynamics of (1.1 1 is governed by an effective Brownian motion 
with a nonnegative covariance matrix, the effective diffusion tensor or effective diffusiv- 
ity. The calculation of the effective diffusivity, in the case where the velocity field is a 
smooth periodic field, requires the solution of an appropriate boundary value problem, 
the cell problem, together with the calculation of an integral over the period of the velocity 
field HSl Ch. 13]. Similar results hold for random, time dependent velocity fields [3]. 

Various properties of the effective diffusivity have been investigated. In particular it has 
been shown that the effective diffusivity is always enhanced over bare molecular diffusivity 
for incompressible fiows |20ll5l9j while it is always depleted for potential flows, [36], (see 
[251 Ch. 13] for a discussion). Furthermore, the scaling of the effective diffusivity with 
respect to the bare molecular diffusion a, in particular in the physically interesting regime 
a <^ 1, has been studied extensively in the literature. It has been shown that the scaling of 
the diffusion coefficient with a depends crucially on the streamhne topology |9|15|35|3^ . 
For example, for steady flows with closed streamlines the effective diffusivity scales like 
a, as a ^ 0, whereas for flows with open streamlines (shear flows) it scales like l/u^ [11] . 
Note that for w = the diffusivity scales like cr^. 

There are various physical applications where modeling the noise in equation ( |1.1| as 
delta correlated in time is inadequate. As an example we mention the problem of transport 
of passive scalars in the ocean; in this case the noise comes from the unresolved velocity 



scales which are correlated in time [1]. A simple variant of equation ( 1.1 ) where the noise 
process has a non trivial correlation structure is 

X — v(x,t) + arj (1-2) 

where r] is an Ornstein-Uhlenbeck process with exponential correlation function. 

We will refer to this problem as the coloured noise problem for passive tracers. It is still 
possible to show, using multiscale/homogenization techniques, that the long time/large 



scale behaviour solutions to equation (1.2 1 is governed by an effective Brownian motion 
with effective diffusivity D [T^]. 

On the other hand, there exist various applications where inertial effects cannot be 
ignored. As examples we mention rain formation [3116] and suspensions of biological or- 
ganisms in the ocean. Recently there has been a burst of activity on the theoretical and 
numerical study of various mathematical models for the motion of inertial particles in 
laminar and turbulent flows [1'33'34J . The starting point for many theoretical investiga- 
tions concerning inertial particles is Stokes' law which says that the force F^it) exerted 
by the fluid on the particle is proportional to the difference between the background fluid 
velocity and the particle velocity: 
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Fs{t) cx v{x{t),t) - x{t). 



(1.3) 



Various extensions of this basic model have been considered in the literature, in particular 
by Maxey and collaborators |16ll7ll8ll9l29l37j . 



The equation of motion for a particle subject to the force ( 1.3 1 and molecular diffusion 

is 

Tx = v{x,t) ~ X + aW (1.4) 

where x G M'^, is a c?-dimensional Brownian motion and t is the Stokes number. We 
will refer to this as the inertial particles model. It is possible to show, with the use of 
multiscale/homogenization techniques, that for either steady periodic or time dependent 



random velocity fields, the long time, large scale dynamics of solutions to ( 1.4 1 is governed 
by an effective Brownian motion |10I24I26] . For the calculation of the effective diffusion 
tensor the solution of a boundary value problem is required, together with the calculation 
of an integral over the phase space T'' x M''. 

The effective diffusion coefficient for inertial particles depends in a complicated, highly 
nonlinear way on the parameters of the problems, such as the Stokes number, the strength 
of the noise and the velocity correlation time. Very little is known analytically concerning 
the dependence of the effective diffusivity on these parameters, in contrast to the passive 
tracers case where this problem has been studied extensively. For example, numerical 
simulations presented in |,24i,26j suggest that for incompressible flows the effective diffu- 
sivity for inertial particles is always greater than that of passive tracers, but no proof of 
this result exists. Furthermore, the presence of inertia gives rise to various asymptotic 
limits of physical interest in addition to the relevant distinguished limits for the passive 
tracers model, such as taking a while keeping the Stokes number r fixed, as well as 
the limit where (cr, r) — > simultaneously. 

Since very little is known analytically for the effective diffusivity, one has to resort 
to numerical (Monte Carlo) simulations for its calculation. For the accurate numerical 
calculation of the diffusion coefficient it is necessary to use numerical methods that can 
integrate accurately the stochastic equations of motion over long time intervals (i.e. until 
the system reaches the asymptotic regime described by an effective Brownian motion); 
furthermore, it is desirable that the numerical method is robust with respect to variations 
in the parameters of the problems such as the Stokes number, the strength of molecular 
diffusion etc. In particular, we want numerical methods which perform well when the 
parameters in the equations of motion become either very large or very small. 

The purpose of this paper is twofold: 

i) We propose new numerical integrators for passive tracers and inertial particles, by 
constructing stochastic generalizations of the geometric integrator proposed in |28j . 
The idea is to construct integrators by composing a geometric integrator with the 
explicit solution of a Gaussian stochastic differential equation. This idea is also used 
in molecular dynamics: see \22\ and the references within. 

ii) Having shown the efficiency of the resulting method for the calculation of the ef- 
fective diffusivity for passive tracers, we investigate numerically various asymptotic 
limits of physical interest, including inertial particles and the coloured noise problem 
for passive tracers. In addition, we generalize the invariant manifold result for iner- 



tial particles in cellular flows [29 to the stochastic case (1.4 1, by using stochastic 
averaging. 
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The rest of the paper is organized as follows. In Section 2 we describe the stochastic 
geometric integrator for both the passive tracers and inertial particles cases. Section 3 
contains some theoretical analysis of the numerical method, more specifically the proof of 
convergence of the method and the behaviour of the method in the case of small inertia. 
In Section 4 we summarize various results concerning equations ( |1.1[ ) and ( 1 .4 1 in the 



asymptotic limits cr ^ and t — *■ 0. Finally Sections 5,6 and 7 contain various numerical 
investigations of the relevant asymptotic limits highlighted in Section 4. 

2. Stochastic Geometric Integrators 

In this section we describe the stochastic splitting method for both passive tracers 
and inertial particles. Before doing this we describe a special feature of the Taylor-Green 
velocity velocity field and then describe the numerical method in a general framework 
abstracting this case. 

The Taylor-Green velocity field v — V^^'tgQ ^tg = sina;i sinx2 can be written as 

( — cosxTsinxi \ 
(2.1) 
sinx2Cosa;i / 

Using the product formula sin a cos f3 = ^ sin (a + /3) + ^ sin {a — f3), inline the vector 
field V can be split as follows: 

v{x) = dig{{ei,x)) + d2g((e2, x)), 

where by (•, •) we denote the usual inner product on M^, g{x) — sin a:, and 

The key property of this velocity field is that the vectors dj, Cj are orthogonal for each j: 

(d,,e,) =0, j = I,--- ,n. (2.3) 

With this in mind we proceed to the analysis of the splitting method in the case where 
the velocity field can be written as 

n 

v{x) = J2djVj{{ej,x)), (2.4) 



assuming (2.3) holds. Note that each vector field djVj{{ej, •)) is itself incompressible and 

onsider 

= djVj{{ej,x')), (2.5) 
(e„x^) = 0. (2.6) 



integrable because of (2.3 1. More precisely, if we consider the ODE 

d 
di 



is easy to check that, by (2.3 1, 

d 
dt 

1 here V-^ stands for V-^ = (- ^ , ^) 
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Thus 



(2.7) 



Note also that in the Taylor-Green case the vectors are 2-dimensional but the method 
we describe works for vectors of arbitary finite dimension d. The idea of exploiting the 
splitting ( |2.4[ ) to construct volume-preserving integrators for with cr = is intro- 
duced in 



2.1. Passive tracers 



We now describe the stochastic splitting method for velocity fields with the properties 



(2. 3), (2. 4) in the passive tracers case. We want to solve the equation (1.1 1. Consider the 



flow (j)j{x,t) generated by (2.7 1: 



(pj {x,t) — X + tdj Vj ( {cj ,x)). 



(2.8) 



Note that (j)j{x,t) is area preserving. A numerical approximation of the deterministic 
part of ( [TT) is given by [55]: 



Xk+l ^cj){xk,/^t) 

(j){x,t) = ((/i„ o • • • o (l)i){x,t) 



(2.9) 



As the composition of volume-preserving maps, (j) is itself a volume-preserving map j3S]. 

(2.10) 



To incorporate the stochastic part of equation ( 1 . 1 1 we simply set 



Xk+i = (l){xk,At) + aV Atjk, 
where jk are i.i.d vectors with 71 ~ A/'(0, Id)- If we define the random map 



then 



■ip{x, t, ^) = 4>{x, t) + aVi^, 

Xk+l = ■l}){xk,A.t,"fk). 



(2.11) 



2.2. Inertial Particles 



We now describe the stochastic splitting method for velocity fields with the properties 



(2.3), (2.4 1 in the inertial particle case (1.4 1. A generalization of the strategy from the 
previous subsection is as follows. The inertial particles system (1.4) can be written as a 
first order system 

z = F{z) + YW 

where z — {x,y), and 



/ 



F(z) 



--y 



]djVj{{ej,x)) 



-.W 
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The most straightforward sphtting would appear to be found by writing F{z) 
with 

/ ^ 

ny/r- 
1 , 1 



--y 



= djVj{{ej,x)) y 



and 



Fn+i{z) — a 







W 



which corresponds to adding a Brownian motion in the last step. However, the resulting 
splitting methods lead to restrictions on At/r, something we wish to avoid. To this end 
we set 

/ 1 



and 



1 



Fn+liz) 



{n + 1)t 

where the last step now corresponds to an Ornstein Uhlenbeck process with timescale r 
in y. The deterministic subequations corresponding to this splitting are 

(n + 1)tx-' = djVj({ej, x-')) — x-' 

and the stochastic part is 

(n + 1)tx = — i + aW. 



(2.12) 
(2.13) 



If we take the inner product with Cj in (2.121 and use (2.4) then we obtain 
{n+l)T{ej,x^) + {ej,x^) =0, j = ,n. 

Thus 



(cj, a;-* (i)) = a — b{n + l)Texp 



where 



(n+ l)r_ 

a={ej,x^O)) + {n+l)T{ej,x\0)), b^{ej,x\0))^ 



Note that now equation (2.12) becomes 

{n + 1)tF + = fj{t), j = 1, • • • , n 
which can be solved explicitly up to quadratures to give the result 



y'(^) =2/^(0) exp 



(71+1)^ Jo 
t 



y^{s)ds, 



(2.14) 



(n + 1)t 



exp 



jt-sy 

{n + l)T 



fj{s)ds. 



where 



fj{t) = djVj{{ej,x^{t))) = djVj a - b{n + l)rexp 



i^^+l)T 



(2.15) 



We denote by (f>j{x,y,t) the flow generated by (2. 14), (2. 151 to obtain a first order 
integrator for the noise free dynamics: 

(l>{x, y, At) = o • • • o cj)i){x, y, At). 



We now take into consideration the stochastic part of ( |1.4| . From (2.131 we have 

^ (2.16a) 



(n+ 1)7? 
-1 

-y 



y, 



(n+l)T- 

It is possible to solve this SDE explicitly to give: 

x{t) 



y{t) 

where A(a:, y, t) is defined by 
Kx,y,t) = 



^Xo<l>{x{0),y{0),t) + g{^,j,t) 



x + [ 1 — exp 



V 



yexp 



(n+ 1)t 

-t 



(n+ 1)t 



and g{'y,C,t) describes the noise. It is given by 



gi^n.t) 



where 7, ^ are i.i.d vectors with ^ ~ A/'(0, /) andp^ 



a^ + d^ = (7^ [t- 2(?i + l)r(l - e" 



(n + l)r 



(1-e- 



/3a = 
/3 



a^^?{n+l) 



1 e 



2 _ {n+l)a^ 



The split-step approximation for incrtial particles is given by 

^{xk,yk,At,£„j), 



Xk + l 

Vk+i 



■^P{x, y, 7) = A ° <i>{x, y, t) + g{^, 7, t) 



where 

^ For details on how to calculate a, (3, 5 we refer to [27 
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(2.16b) 



(2.17) 



(2.18) 

(2.19) 
(2.20) 



Before we proceed to the next section we briefly discuss how we solve (2.141. The integrals 
in (2.141 can be calculated with some high level quadrature. However we take advantage 
of the fact that we can substitute for (t) in the (t) equation and convert the double 
integral into a single integral to obtain 



{t) = x' (0) + Vrj/J (0) (^1 - exp 
^ {t-s) 







1 — exp 
J/-' (i) =2/^(0) exp 



(n+ l)r 
t 



(n + 1)t 



(n+ l)r 
fj{s)ds, 



exp 



(2.21) 



(n + 1)t 



fj{s)ds. 



In order to calculate the convolution integral arising in both cases we make the substi- 
tution q = e to obtain 



exp 



(n+ 1)t 



fj{s)ds = (n+ l)Te 



djVj I a 



b{n + l)r 



dq 



Now we approximate the integral on the right hand side by a simple Euler method to 
obtain 



e ("+!)'■ 
where 



djVj I a - 



b{n+ 1)t 



w (1 - e ("+!)- )dji;j(a - br^n + l)). 



(l_e ('^+^>-)djVj{a-bT{n+l)) = {l-e <"+i)- )djWj((ej, x^'(O))) 

while the other integral in the x equation in (2.21 1 is approximated using the trapezoid 
rule. Thus we obtain the following approximation to (f)j(x,y,t), denoted by (f)j{x,y,t), 
where 



and 



with 



Hj{x,y,t) 



<Pj {x, y,t) ^ X{x, y, t) + fij{x, y, t) 

(n + 1)t(1 - e^T^^)djVj{{ej,x)) + hj{x,y,t) 
(n+ 1)7^(1 - e-T^)d,Vj{{ej,x)) 



hj{x,y,t) ^ dj- Vj{{ej,x))+Vj{{ej,x+{l-e ("+1'-)?/)) 
We find this effective in practice, especially for small r. 

3. Analysis of the numerical method 

3.1. Convergence of the stochastic splitting method 

In this subsection we present a result regarding the strong order of convergence for the 
stochastic splitting methods. Note that we can write (1.1 1 and (1.4 1 as a system of first 
order SDE's namely: 

F{z) + Y.W, (3.1) 
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where z e W and I = d or I = 2d, F{z) = Y.'lF,{z), and S e E'^'' and W is an 
d-dimensional standard Brownian motion . We now state a theorem concerning the con- 
vergence of the numerical method. 

Theorem 1 /i^/ Let a;„ he the numerical approximation of ( |3.1| for the stochastic split- 
ting method at time nAt, where F{x) — '^^^iFi{x), Fj E C'^(M', M'). Suppose that 

E{\xo\^)<oo, (3.2a) 
E(|xo-2/o^f )i/'<ifiAii/2^ (3.2b) 
\F,ix)-F,{y)\<K2\x-y\, J = !,■■■ ,n (3.2c) 

where the constant Ki^K^ do not depend on At. Then 

\ 1/2 

E sup |a;(fcAi) - Xk\^ < K2{T)At, V T > 0. (3.3) 

a<kAt<T J 

Note that the strong order of convergence is 1, as in the case of the Euler-Marayama 
method with additive noise il3i. 



3.2. Splitting method in the case of small inertia 

In this subsection we investigate the behaviour of the stochastic sphtting method in 
the case of small inertia. Our interest is in studying the behaviour of the method as 
T 0. When we send r to 0, while keeping At fixed, we recover the solution of the split- 
ting method for the passive tracers problem. This is the content of the following theorem. 



Theorem 2 Let Xi^t(T), x'^^{T) denote the numerical approximations to equations (1.1 ), 



(1.4) obtained at time T — nAt using the stochastic splitting method. Let the assumptions 



of Theorem^ hold. Then: 



(E(||xA.(r) - xl,{T)\\')f' < C{T)^ (3.4) 



Proof: We need to introduce the operators and Py 





(x\ 




(x\ 






= Py 













y, 



where x,y E M'*. In Section 3 we defined the maps ?A(a;, y, At, 7, ^), -0(3:, At, 7) and ex- 
pressed the numerical solution for both passive tracers and inertial particles through 

passive tracers, Xk+i — 'ip{xk,At,jk) 

inertial particles, {xk+i,yk+i) = -ipixk^yk, At,^k,7k)- 

We now set Ck = Xk — Xk to obtain 
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efc+1 ^ Pxi^{xk,yk,At,jk,^k) - -(/"(xfe, Ai,7fc) 

= Pa;(A o (j)){xk,yk,/^t) + Pj;g{^knk, At) - (j){xk,At) - a^/Atjk 

At 
(n+ l)r 



= Px(t>{xk,yk,At) + ^rPy 1 - cxp ( -y— ^^^^ ) (l){xk,yk,At) ~ (j){xk,At) 

+ Px9{£.k , Ik , At) - aVAijk ■ 
We take norms and use the triangle inequality to obtain 

||efc+i|| < \\Px^{xk,yk,At,) ^ (l){xk,At)\\ + y^\\Py^{xk,yk,At)\\ 
+ \\Px9khk,^k,At) ~ ay^t^kW- 
Using Lemma [12] from the appendix we obtain 

||efc+i|| < {l + KAt)\\ek\\ + CV^+\\Pxgkhk,^k,At)-a^tjk\l (3.5) 

and if we take expectations and use Lemma[T3]from the appendix together with Jensen's 
inequality for the noisy part we conclude 

E(||efc+i||) < {l + KAt)E{\\ek\\) + {M + C)V^. 

We now use the discrete Gronwall inequality and set Ci — M + C to obtain 



{l + KAt)^-l 



since (1 + LAf)" < e"-^'^*. Hence we deduce that 

kKAt 



KT 



1 



KAt ' 



efcll) <Cie« 



At' 



(3.6) 



We now use (3.5 1 again by taking squares and then expectations to obtain 

E(||efc+i|n < (1 + XAt)2E(||efc|n + 

+ 4 [Cx/7 + E(||P,.g(6,7fe, Ai) - aAt-fk\\)] E(||efc||) 
+ 2CV^E{\\Pxg{^k,lk,At) - <jAt-fk\\) 
+ E{\\Pxg{^k,lk,At) - aAt-fk\f). 

where in the second line we have used the fact that KAt < 1. We can now use the 
equation (3.6 1 together with Lemma 13 from the appendix to obtain 

r 



E(||efe+i||^) < {I + LAt)E{\\ek\n + Mt + Cie 
By applying the discrete Gronwall inequality we conclude that 



KT 



{E{\\xAt{T)-xUT)r)f'<C{T) 



At 



At 



□ 



4. Relevant Asymptotics Limits 



In this Section we describe the asymptotic limits which guide our numerical experi- 
ments. We start by presenting results concerning the vanishing molecular diffusion limit 
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of the Taylor-Green velocity field and the shear flow in the passive tracers model (1.1 1 
. We also present a result regarding the effective diffusive behaviour of passive tracers 
driven by colored noise. 

We then present results for inertial particles, giving a new bound for the effective 
diffusivity of inertial particles for the shear flow. The case of small inertia is also studied 
and a result is then presented for the relation between the effective diffusivity of inertial 
particles for small inertia and the effective diffusivity of passive tracers. Finally, in the 
last subsection we present a modified passive tracers model obtained from averaging the 
inertial particles model in the case a = ^fr and r ^ 1. Throughout this section we use 
the following definition for the effective diffusivity. 

Definition 3 The effective diffusivity matrix is defined (when it exists) as 

^ ^ {{x{t) - xm ® {x{t) ~ x{m 

t^oo 2t 



where x(t) is the solution of the equations of motion (i.e. (1.1 1, (1.2 I or (1.4)j and (•) 
denotes ensemble average. 

It is possible to prove rigorously for all the problems that we will consider in this paper. 



namely equations (1.1 1 (1.2), (1.4 1 with time independent, periodic and incompressible 
velocity fields, that the effective diffusivity exists. More precisely, the rescaled process 

x'{t) := ex{t/e^), 

where x{t) is the solution of the equations of motion, converges weakly (as a probabil- 
ity measure over the space of continuous functions) to a Brownian motion W{t) with 
covariance matrix /C: 

x\t) ^ V2]CW{t). (4.2) 



However we use (4.1 1 to compute /C, by Monte Carlo techniques. 



4.1. Passive tracers 



In this subsection we present results regarding passive tracers in the small molecular 
diffusion limit together with a result concerning passive tracers driven by colored noise. 

4.1.1. The small molecular diffusion limit 

In this subsection we describe two results concerning the behaviour of the effective 
diffusivity in the small molecular diffusion limit for the passive tracers. We study both 
the Taylor-Green velocity field and the shear flow. 

Taylor- Green Velocity Field 

Result 4 \819\15^ Let x{t) he the solution of the passive tracers with v{x) given by equa- 



tion (2.1 1. Then the following results holds in the case a 
where I2 is the two dimensional unit matrix. 

Note that in the case where v(x) — then the effective diffusivity matrix is 
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2 

Thus the relative enhancement of the effective diffusivity for the Taylor-Green field over 
the bare molecular diffusivity is unbounded in the limit cr ^ 0. 



The convergence result (4.2) does not provide us with any information concerning the 
relevant time scales, in particular the time needed for the systems to reach the asymptotic 
regime which can be described through an effective Brownian motion. The scaling of the 
diffusive time tdiff with respect to the molecular diffusivity a was studied in [5^ for 
various types of incompressible flows, for the passive tracers problem. 



Result 5 Let x{t) be the solution of passive tracers with v{x) given by equation (2.11 
and tdiff be the time it takes for the particle to start behaving diffusively. Then 

tdiff ~ < 1- 

Shear Flow 

Another flow of interest is the shear flow 

v{x) = f ° I • (4.3) 
\ sinxi / 

For this flow explicit calculation gives the following theorem: 

Result 6 fWUM^jj Let x(t) be the solution of passive tracers with v{x) given by equation 



(4.3). Then the effective diffusivity matrix is given by 

2 






JC{a) = I ^ ^2 1 I . (4.4) 

Note that the result shows that the effective diffusivity is unbounded, in absolute terms, 
in the second component of the system as cr ^ 0. This remarkable effect arises from 
ballistic transport over long distances, slowly modulated by molecular diffusion. 

4.1.2. Passive Tracers Driven by Colored Noise 

In various applications it is sometimes of interest to consider passive tracers driven by 
coloured noise (see U and references within). The equations of motion are then 

(TTI 

x^v{x) + -^, (4.5a) 
r,= -'l + ^W. (4.5b) 



Techniques from homogenization show that equation (1.1 1 is recovered in the limit ^ ^ 



It can be proved then that under appropriate assumptions on the velocity field the 



effective behaviour of x governed by (4.5) is Brownian motion with an effective diffu- 



sivity matrix ]C{a,S). For the derivation of the effective equation using the backward 
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Kolmogorov equation for such systems we refer to |12j while the same derivation from 
the point of view of the Fokker -Planck equation can be found in [¥ . The following result 
describes the limit 6^0. 



Result 7 1121 Let JC{a,S) be the effective diffusivity matrix for x(t) governed by (4.51 
and K,{a) the effective diffusivity matrix for x(t) governed by (1.1 1. Then, for a fixed and 
(5< 1, 

IC{a,5)^ICicr) + C{a)5 + 0{S^/^). (4.6) 



4.2. Inertial Particles 



In this subsection we present some relevant results for inertial particles. 



4.2.1. Shear Flow 

We can use Result [6] to prove the following result for inertial particles: 
Result 8 /U/ Let x{t) be the solution of (1.4 1 with v{x) given by (4.3 1. Then the effec- 
tive diffusivity matrix is given by 



JC{a,r)= \ 2 

/C22(o-, r) 



where, for some Ci,C2 independent of a and t, 

\IC22{cr,T) - /C22(a-)| < 



2Ci 



C2T. 



(4.7) 



Here /C22(o') is given by (4.4 1 as 



/C22(o-) = 



Note that equation (4.7 1 shows that for a <C 1 and r = 0{1) , /C22(o', t) behaves like 
1^22(0') since in this case /C22(ct) grows like l/cr^. This explains the numerical results in 
[24j . where it was shown that K. is effectively independent of r in the inertial particles 
case. 



4.2.2. The case of small inertia 

In this subsection we study the case of small inertia. We have the following theorem, 
using techniques from |^ : 



Theorem 9 ^121 Let x, x'^ solve the stochastic differential equations 

X — v{x) + crW, 
TxT — v{x'^) — + aW. 

where v G C{T'^,W^). Then 



E sup ||x^(i) -x{t)\[ 

0<t<T 



1/2 



< i^Wrlog ( ^ +2 )e 



LT 



(4.8a) 
(4.8b) 
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Theorem [9] gives us only pathwise information, but it does not reveal the relation between 
the effective diffusivity of inertial particles and passive tracers in the small inertia limit. 
The following result relates the effective diffusivity of passive tracers with the one of 
inertial particles in the small r regime 



and K,{ct) the effective diffusivity matrix for x(t) governed by (1.1 1. Furthermore assume 
that V • u = 0. Then for a fixed and r ^ 1, 



Result 10 Ii^6] / Let /C(cr, t) be the effective diffusivity matrix for x{t) governed by (1.4 1 

ty matrix for x(t) governt 
d and r ^ 1, 



4.2.3. The case a — ^/t 

In this subsection we present a result for the case where a — yfr and t ^ using 
formal asymptotics arguments. It is exactly this relationship between tr, r that makes 
such a treatment possible, since then the leading order operator is ergodic and thus we 
can apply the techniques of stochastic averaging. 

The next result follows from formal perturbation arguments, similar to those used in 
Part II of [33]; details are given in the appendix. 

Result 11 For r 1 and t — 0(1) the dynamics of the inertial particles model 

Tx^ = v{x) - + y^rW (4.9) 

are approximated by the modified passive tracers model 

X ^ v{x) - t{Vv{x))v{x) + y/rW. (4.10) 

In the case r = the modified passive tracers model is simply the equation for Lan- 
grangian trajectories, while in the absence of Brownian motion it is precisely the two-term 
approximation of the invariant manifold found in |29j . 



5. Numerical Investigations: The Vanishing Molecular Diffusion Limit 

In this Section we investigate the performance of the stochastic splitting method in 
the vanishing molecular diffusion limit. We study its behaviour for both passive tracers 
and inertial particles; and we use both the Taylor-Green velocity field and the shear flow. 
The objectives of our investigation in this section are as follows: 

i) To compare the stochastic splitting method with the Euler method in the passive 
tracers and the inertial particles case. 

ii) To use modified equations to show that the Euler method is not suitable for the 
passive tracers case. 

iii) To apply the stochastic splitting method in a slightly different setting, the case of 
passive tracers driven by colored noise, thereby obtaining new information about 
this problem. 

IV ) To verify Result |8] concerning the behaviour of inertial particles under the shear flow 

in the vanishing molecular diffusion limit, 
v) To obtain new results for the vanishing molecular diffusion limit in the case of inertial 
particles. 
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5.1. Passive Tracers 



In this subsection we study the behaviour of the stochastic sphtting method in the 
vanishing molecular diffusion limit for passive tracers. Note that the behaviour of the 
effective diffusivity is analytically known for the shear flow (Result |6| and the Taylor- 
Green velocity field (Result |4]). This provides us with a good testing ground for the 
stochastic splitting method. Note that in the case of shear flow the Euler method and 
the stochastic splitting method are the same and thus we do not present any results for 
the shear flow in the case of passive tracers. 

We now present some numerical results concerning the passive tracers problem as 
cr ^ for the Taylor-Green velocity field. In this case the effective diffusivity matrix is 
diagonal with the diagonal elements scaling like a. 

In Figure [T| we plot the phase plane of (1.1 1 for a — 10^^ with the use of the two 
different numerical methods. The realization of the noise is the same for the two methods 
and we have integrated up to time T = 10** with timestep At = 10~^. 




3^ 



(a) Euler method (b) Stoehastic splitting method 

Fig. 1. Phase plane for the two different methods. 

It is clear that the behaviour of the particle is drastically different. In the case of the 
Euler-Maryama method the particle appears to be much more diffusive than in the case 
of the stochastic splitting method. 

We now compare the two methods in the case of zero noise, since that will help us to 
understand the different behaviour in the small noise regime. In Figure [2] we draw the 
phase plane in the absence of noise. 

In the absence of noise the system is Hamiltonian and thus solutions to ( |1.1| should 
follow closed orbits. The stochastic splitting method maintains this property of the system 
, being a volume preserving method, as shown in Figure [2]3. On the other hand the Euler 
scheme fails to maintain this property of the system, and instead solutions spiral out |llj 
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(a) Euler method 



(b) Stochastic splitting method 



Fig. 2. Phase plane for the two different methods in the absence of noise. 

- see Figure [2^. When a is small the spiraling effect in the Euler method leads to faster 
escapes from the cell than the stochastic splitting method and hence to overestimated 
effective diffusivities. 

In Figure[3]we plot /Cn as a function of the molecular diffusivity cr for the two different 
methods. For both of the methods we have used N = 10'^ realizations integrating up to 
T = 10^ with time step At ~ 10~^. We can clearly see that for large values of molecular 
diffusivity the two methods agree but as a gets smaller the Euler method fails to capture 
the results predicted by theory. The stochastic splitting method, however, agrees with the 
theory. If we fit our data for values of cr G [0.005, 0.1] we find that the effective diffusivity 
grows like c*a"-, where a — 1.0946, c* — 1.0575, while the theory predicts that a ~ 1.0. 

We now reduce the time step to At = 0.005 to obtain Figure |4] and consider only 
the stochastic splitting method. We can clearly see that the numerical results agree 
almost perfectly with the the theory for the region of cr G [0.005,0.5]. If we now fit 
the data for cr G [0.005,0.1] we find that the effective diffusivity grows like c*cr°, where 
a = 1.0579, c* = 0.9269, which agrees almost perfectly with the theory. 

For values of cr smaller than 5-10^^ we see that the effective diffusivity does not behave 
as the theory predicts, even for the splitting method. However this should not come as 
a surprise, since as we see in Result [5] the time for the particle to be diffusive is of order 
1/cr^ and thus integrating up to T = 10^ is not enough since the particle has not yet 
reached its diffusive regime. 

Before rejecting the Euler-Maryama method for our problem we will make one last 
comparison. We expect that the volume preserving method is roughly 3 times slower 
than the Euler method since it involves 2 extra steps in order to compute the value of 
the solution at each timestep. Thus, a fair comparison is to compare the Euler method 
over the splitting method when the former uses a timestep that is one third of the 
latter. So, we use a timestep of At — 10"^ for the volume preserving method and then 
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Fig. 3. Effective diffusivity as a funetion of a for the two different methods. 



Fig. 4. Effective diffusivity for cr <C 1 using the stochastic splitting method. 

calculate the effective diffusivity for a = 10^^ using N = 10"^ trajectories. We repeat the 
calculation with the Euler method. We then decrease the timestep by successive factors 
of two until our final timestep becomes At = 10~'^/2^. The results are plotted in Figure 

H 

The failure of the Euler method is once again obvious. Even when we use the a timestep 
16 times smaller than the one used in the the splitting method, the calculated effective 
diffusivity is 5 times larger than the correct value, which we calculate with the splitting 
method. 

The Euler method not only forces the particle to behave more diffusively but also to 
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time- step 

Fig. 5. Comparison of the two methods for u = 0.01 and different timesteps. 

reach its diffusive regime much faster than the theory predicts. This is clearly exhibited 
in Figure [6] where we plot the effective diffusivity as a function of time. We know from 
Result [s] that the time it takes for the particle to start behaving diffusively is 0(l/cr^). 
This is indicated by the vertical line in Figure [6] The first line from the top corresponds 
to the Euler-Maryama method for At — 10~^/4 and the two below for At — 10~^/8, 
At = 10~^/16, while the last line corresponds to the stochastic splitting method for 
At = 10^^. It is clear that the particle reaches its diffusive regime faster than it should 
when the Euler-Maryama method is used. 

5.2. Modified equations 

In this subsection we study the Euler method with the use of modified equations. 
Modified equations is a widely used method for backward error analysis |11| for ordinary 
differential equations. In the case of stochastic differential equations the derivation of 
modified equations is more complicated, because the semigroup governing expectation 
propagation is not invertible. However, some limited work is available and we employ a 
form of weak backward error analysis from [30j . 

Consider the following SDE 

dX = f{X)dt + a{X)dW{t), X(0) = Y (5.1) 

where X G M'^ and W a d-dimensional standard Brownian motion. Let us consider a 
numerical approximation Xq, Xi, ■ ■ ■ Xn of weak order 1 

|E0(X„) - E(t){X{nAt))\ = 0{At), < nAt < T, (5.2) 
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Fig. 6. 



<(:ri(t)-a:i(0))^) 



as a function of time for the two different methods. 



for in a space of smooth test functions [13 . We would like to modify the SDE (5.1 1 to 
define a process X that better describes the numerical approximation Xn in the sense 
that 

|E^(X„) -E(/)(X(nAi))| = ©(At^), < nAt < T. 
We define X as the solution to the modified SDE 



dX 



f{X + f{X))At 



dt 



(T{X) + d{X)M dW{t), X{Q)^Y 



(5.3) 



where /, a are smooth functions to be determined. In [30 /, a were derived in the case 
of the Euler method with additive noise. For our equation (4.8 1) the modified equation 
is 



v{x) 



At 



{Vv{x))v{x) 



Av{x) ] dt + a [1 



At. 



Vv^{x) W. (5.4) 



Note that the correction proportional to a^At in the drift is related to the presence of 
noise in the problem since in the absence of noise the modified equation for the Euler 
method would only contain the (\7v{x))v(x) correction |llj . 

If w = V^^* we apply Ito's formula to the stream function ^'(x) where x satisfies (1.1 1 
and we find 



-A^' + M. 



Here the integral of M is a zero mean martingale. 

In the case of the Taylor-Green fiow A^* = — 2^*, so equation (5.5 1 becomes 



M. 



(5.5) 



(5.6) 



Thus the mean value of the stream function decays like e"^!" . We now apply Ito's formula 
to the stream function V'(^) for the Taylor-Green fiow for x satisfying the modified 
equation (5.4 1. We find that 
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At, 



2 (cos^ xi + cos^ 3:2)^' — (T^^'(l + Atcosxi cosa;2) 
H , — (cos^ x-Y cos^ X2'^ ~ ^' ) + Ma*, 



(5.7) 



where the integral of Ma* is again a mean zero martingale. 



Note that in the case At = equation (5.7 1 becomes (5.6). Furthermore the first term 
on the right hand side of (5.7l does not depend on a. This term when Ai ^ C'(o'^) causes 



the spiraling effect seen in Figure [2^. Furthermore, it drastically changes the behaviour 
of the mean value of the Hamiltonian ^' as a function of time as seen in Figure [7] where 
the mean value of the Hamiltonian as a function of time is plotted for the two different 
methods, together with the real solution. As we can see, in the case of the Euler method 
the mean Hamiltonian decays much faster than the theory predicts, while in the case of 
the stochastic splitting method it decays at the right rate. 




Fig. 7. Mean value of the Hamiltonian as a function of time, for At = 10 ^, = 10 ^. 



5.3. Passive Tracers Driven by Coloured Noise 

In this section we study numerically the motion of passive tracers in the Taylor- Green 
velocity field subject to coloured noise using a slightly altered stochastic splitting method. 
More precisely, the deterministic steps are the same as in the case of passive tracers driven 
by white noise. However, when we add the noise we add an exactly sampled OU process, 
using similar arguments as in the case of small r for inertial particles to sample exactly 



the OU process and its integral in equation (4.5 1, instead of just adding white noise 
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(a) 5 = 10-1 (b) <5 = 1 

Fig. 8. Effective DifFusivity as a function of <j for the two different metfiods. 



In Figure [s] we plot the effective diffusivity as a function of cr for 5 = 10~^, and S = 1 
for the two different methods. In our calculations we have used a time step At = 10~^, 
with final integration time T = 10^ for N = 10'^ paths. 

In the case of 5 = 10~^ we expect that the effective diffusivity for the coloured problem 
should be close to the one of the white problem (Result |7]) and thus it should go to zero 
as fj ^ 0. As we can see this is the case for the stochastic splitting method, but not 
for the Euler method. This should not be a surprise since as we have already seen in 
the case of white noise the Euler method does not capture the right behaviour for the 
effective diffusivity. We also see that similar behaviour is seen for S = 1. In particular, 
the numerical experiments suggest that the effective diffusivity when cr ^ 1 is essentially 
independent of S. 



5.4. Inertial Particles 

In this subsection we study the performance of the stochastic splitting method for the 
case of inertial particles,, and we also study the dependence of the effective diffusivity 
on the various parameters of the problem. We study both the shear flow and the Taylor- 
Green velocity field. In doing this we should keep in mind that the behaviour of the 
effective diffusivity is known for the shear fiow (Result |8| but not for the Taylor-Green 
velocity field. 
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5.4.1. Shear Flow 

In this section we present some numerical results concerning the inertial particles prob- 
lem in the vanishing molecular diffusion limit for the shear flow. Note that, in contrast 
to the passive tracers, the Euler and the stochastic splitting method differ in the inertial 
case. In Result |8] we have shown that in the molecular diffusion limit the /C22 element 
of the matrix diverges like l/cr^. Thus, this is the scaling that we expect our method to 
capture as ct ^ for r of 0(1). 

We now start our investigation by computing effective diffusivities with the two dif- 
ferent methods. In Figure [9] we compare the effective diffusivities for the two different 
methods. For this figure the time step used was At = 10^^ and we integrated for time 
T = 10^, and for N — 10'^ realizations. If we now try to fit the data for a e [0.01, 0.5] we 



Fig. 9. Effective diffusivity as a function of cr for the two different methods for the shear flow. 

find that the effective diffusivity behaves like ba"^ where b = 1.077, a = —1.9665 for the 
stochastic splitting method and b = 1.0892, a = —1.9706 for the Euler method and both 
of them agree with what is predicted by Result |8] 

5.4.2. Taylor-Green Velocity Field 

In this section we present some numerical results concerning the inertial particles 
problem for the Taylor-Green flow in the limit of small diffusion. 

Note that in this case no analytical result for the behaviour of the effective diffusivity 
is known, in contrast with the passive tracers case. However in [26' numerical evidence 
was presented indicating that the presence of inertia enhances further the diffusivity. 
Thus we expect that if the effective diffusivity goes to zero in the a — > limit, this would 
happen no faster than linearly, since this is the case for the passive tracers. 
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(a) Euler method (b) Stochastic splitting method 

Fig. 10. Phase plane for the two different methods. 



We start our investigation as we did in the case of passive tracers by comparing the 
two methods pathwise. We choose the value of a — 10~^ and integrate for T — 10"* with 
time step At =10^^ using the same noise reahzation in the y equations, to obtain Figure 

As we see in Figure [TO] the quahtative behaviour of the solutions is the same, unlike 
the passive tracers case. We see in Figure [Tl] that the effective diffusivity is the same for 
the two methods with r = 1, for values of cr e [0.01, 1]. The effective diffusivity has been 
calculated with time step At = 10~^, final integration time T = 10^, and for N = 10^ 
realizations. 

In Figure [12] we compare the effective diffusivity of passive tracers and inertial parti- 
cles in the small diffusion regime, calculated in both cases with the stochastic splitting 
method. The results are consistent with the conjecture in [26 that the effective diffusiv- 
ity of inertial particles is always greater than the passive tracers one. Note also that the 
behaviour of the effective diffusivity as a function of a in the case of inertial particles is 
highly nonlinear, unlike the passive tracers case. 

Concluding this investigation we see that for the case of inertial particles under the 
Taylor-Green velocity field, it seems that there is no advantage in using the stochas- 
tic splitting method over the Euler-Maryama method in the small molecular diffusivity 
regime, for r of 0(1). 
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Fig. 11. Effective diffusivity as a function of a for the two different methods for the Taylor-Green flow, 
T = 1. 



Fig. 12. Comparison of the effective diffusivity: passive tracers vs inertia! particles. 

6. Numerical Investigations: The Small Inertia Case 



In this section we study the effect of small inertia on the problem ( 1 .4 1 . We know from 
Result 10 that the first order correction to the effective diffusivity matrix is 0{-\/t) and 



that is what we would like our numerical method to reproduce. 
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As we have previously mentioned, in this regime the Euler method will fail to give us 
an accurate calculation of the effective diffusivity in any reasonable computing time. The 
reason for this is that we require At = 0{t) in order to avoid numerical instability. This 
makes the use of the Euler method impractical for very small values of r. Thus from now 
on, we will use the stochastic splitting method to investigate the small inertia limit. 

Our objectives in this section are: 

i) Verify Theorem [2] numerically. 

ii) Study the asymptotic behaviour of the splitting method, for r ^ 1. 

iii) Study the behaviour of the effective diffusivities for t ^ 1. 



6.1. Shear Flow 



In this subsection we study the effective diffusivity of inertial particles under the shear 
flow in the small inertia regime. As we have already seen in Section 4.2 the stochastic 
splitting method for fixed At maintains the pathwise convergence to passive tracers as 
T ~> 0. 

This property of the method is illustrated in Figure 13 



where we plot supj^ 



<nAt<T I 



Fig. 13. Log-log plot of the error as a function of r. 

x]-i\\ as a function of r, where At ~ 10^'^, T = 1 and cr = 1. In this case Xn, is the 
numerical approximation for passive tracers and inertial particles moving in the shear 
flow respectively, both calculated using the stochastic splitting method. 

In order for this comparison to be meaningful we have used the same noise realization 
in the x equations for the passive tracers and the inertial particles (see the proof of 
Lemma 13 in the appendix) 
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Fig. 14. Effective Diffusivity as a function of r for tr = 1. 



If we now try to fit the data for r g [10^6, 10^4] in Figure 13 we find that the error 
reduces like where (3 = 0.4955, which is what we expected from Theorem [2] We now 
proceed to investigate what happens for the effective diffusivity of the inertial particles 



in the small inertia regime for the shear flow. In Figure 14 we plot the effective diffusivity 
of the inertial particles for different values of r and cr = 1. For the calculation of the 
effective diffusivity we have used N = realizations with final integration time T — 10^ 
and time step At = 10^^. The results are consistent with Result 10 



6.2. Taylor-Green Velocity Field 



In this subsection we study the effective diffusivity of inertial particles under the Taylor- 
Green velocity field in the small inertia regime. In Figure 



15 



we plot supi<„At<T lla^r, 



as a function of r. Again we have used the same noise m both the passive tracers and 
the inertial particles equation, as in the case of the shear flow. We have used time step 
At = 10~^ flnal integration time T — 1 and a = 1. If we now flt the data of the flgure 
for values of r G [10^^, 10"^] we find that the error reduces like t^, where /3 — 0.5266 
which is close to the prediction from Theorem [2] 

We now proceed with investigating the effective diffusivity in the small r regime for 
the Taylor-Green velocity field. In Figure [I6] we plot the effective diffusivity for different 
values of t for a- = 0.1. We have used a final integration time T — 10* and N = 10^ 
iterations and time-step At = lO"'^. As we can see the effective diffusivity behaves in the 
expected way (Result 10 1, since as r ^ it converges to that of passive tracers. Also, if 



we fit the data for values of r e [0.01,0.5] we find that the effective diffusivity reduces 
like t", where a = 0.504, which again is in agreement with Result \T0\ 
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Fig. 15. Log-log plot of the error as a function of r 




rf 10-' II* 10^ 10-' 11' It" 

T 

Fig. 16. Effective diffusivity as a function of t and At for 
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7. Numerical Investigations: Vanishing Molecular Diffusion and Small 
Inertia 



In this section we study the hmit of both r and a going to zero with the special scahng 
(7 = y^r (see Result 11). In this case we have shown using formal asymptotic arguments 
that the inertial particles for in this particular scaling can be approximated by equation 
( |4lO l: 

X = v(x) — T{yv(xy)v(x) + ^/tPi- 

We will refer to this model as the modified passive tracers model. The objectives of our 
investigations in this Section are: 

i) Extend the stochastic splitting method to study the effective diffusivity for the 
modified passive tracers model. 

ii) Deduce properties of the effective diffusivity in this regime. 

iii) Compare the effective diffusivity of the modified passive tracers model with the 
effective diffusivity of the original passive tracers and of inertial particles and verify 
the validity of the modified passive tracers model. 



7.1. Shear Flow 



In the limit of vanishing molecular diffusion and small inertia the modified passive 
tracers model is trivial since 

{\/v{x))v{x) = 

for the shear flow and so there is no first order invariant manifold correction. Thus the 
modified passive tracers model reduces to 

X — v{x) + ^frW . 

This is precisely the passive tracers models when a = y/r and the behaviour of the 
effective diffusivity is analytically known. 



7.2. Taylor-Green Velocity Field 

In this subsection we study the Taylor-Green velocity field in the limit of vanishing 
molecular diffusion and small inertia. The limit in this case is no longer trivial since 

1 / sin2a-i 
{Vv{x))v{x) = - 

^ \ sin 2x2 

We study the effective diffusivity for this problem numerically. Since we expect the 
Euler method to perform poorly on the modified passive tracers equations, we use a 
generalization of the stochastic splitting method. We perform the first two splitting steps 
as we did in the case of passive tracers. We then add a third deterministic step, where 
we add the correction — r(Vf and finally add the noise. 

Before we proceed with our numerical investigations, we try to get some insight into 
the effect of the correction term on the individual trajectories. One way of doing this is 
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(a) Original passive tracers for cr^ = r (b) Modified passive traecrs model 

Fig. 17. Phase plane for the modified and the original passive tracers model , T = 10'*. 



to study the time derivative of the stream function ^{xi^ X2) — sinxi sin 2:2 for equation 



(4.101. Using Ito's formula we obtain 

^ ^ -T'l!{cosxl+cosxl) -T-^ + M, (7.1) 
where the integral of M is a mean zero martingale. Note that if we compare this equation 



with (5.6) for a — ^fr we see that we have an extra term in the drift, which as we have 
already seen in the analysis of the Euler method, is responsible for a spiraling out effect 
within cells. We thus expect the effective diffusivity of the approximate model to be 
greater than the effective diffusivity of the passive tracers for the same value of the 
molecular diffusion coefficient a. 



In Figure 17 we plot the phase plane for the two different models for final integration 
time T = 10^ with time step Ai = 10~^ and for value of r = 10"'^, using the same 
noise to drive both equations. It is clear that in the case of the approximate model, the 
particle behaves more diffusively than in the passive tracers model. In Figure [18] we plot 



the effective diffusivity as a function of a = y/r for the approximate model (4.10 1 as well 
as for the original inertial particles and passive tracers. 



As we see in Figure 18 the effective diffusivity of the approximate model is larger than 
the effective diffusivity of the original passive tracers model, as expected. Moreover, the 
effective diffusivity of the inertial particles in the limit t — > 0, with a = \/r is greater, 
than the effective diffusivity for both modified and original passive tracers. However, the 
effective diffusivity of the modified passive tracers model captures the full model in the 
case cr^ = r, T ^ 0, which verifies the validity of the approximate model, derived in 
Result [m 
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Fig. 18. Effective Diffusivity as a function of cr = ^Jt for the two different passive tracers models. 

8. Conclusions 

The problem of numerically calculating effective diffusivities in the vanishing molecu- 
lar diffusion limit was studied in this paper. Both passive tracers and inertial particles 
models have been considered. A stochastic splitting method has been proposed that takes 
explicitly into account the features of the equations for passive tracers in the absence 
of molecular diffusion, building on ideas in 28J. Using this method to calculate the ef- 
fective diffusivity we found excellent agreement with the existing theoretical predictions. 
A series of numerical calculations were performed for passive tracers to compare the 
stochastic splitting method and the Euler method; these numerical tests exhibited the 
superior behaviour of the stochastic splitting method. These observations were quantified 
with means of backward error analysis, which revealed the failure of the Euler method 
to capture the essential dynamics for small molecular diffusion. In the case of inertial 
particles for Stokes number r of 0(1), we found no significant advantage of the stochas- 
tic splitting method over the Euler method. However, in the limit of small r we were 
able to prove that for a fixed timestep A< the stochastic splitting method maintains the 
property of pathwise convergence of the inertial particles model to the passive tracers 
one as T ^ 0, for a fixed final time T. This behaviour was manifest in the numerical 
calculation of the effective diffusivities, where the correct limiting behaviour was also 
captured. The combined limit of small inertia and vanishing molecular diffusion was also 
studied in this paper. A modified passive tracers model that approximates the dynamics 
of inertial particles in this regime has been found by means of formal asymptotics. An 
extension of the splitting method was used to calculate the effective diffusivity and the 
results obtained agreed with the theoretical predictions. 
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Appendix A In this appendix we derive the modified passive tracers equation from 
Result 
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We write (4.9 1 as the first order system 



y 

: -v{x) y 



-.W 



(A. la) 
(A. lb) 



The generator C associated with this Markov process is of the form 



T 



where 



Co = v{x) -Vy - y -Vy + - Ay 
Ci^y-Vx 

Note that £o is the generator of an Orstein-Uhlenbeck process in y with mean value 
v{x). The invariant measure of the process is a Gaussian I\f{v{x), \ld) where Id is the d 
-dimensional unit matrix. From now we use {■) p to denote averages with respect to this 
invariant measure noting that C^p — 0. Note that the Fredholm alternative applies: the 
null space of the generator Cq is one-dimensional and consists of constants in x. Moreover 
the equation CqJ = g has a unique (up to constants) solution if and only if 

(.9>p — / 9{y)p{y)dy = o. 



y} denote the solution of (A.ll starting at 



Let X^^y := {x{t),y{t)-x{Q) = x,y{Q) 

{x, y} and let / : M'' x M'' M be a smooth bounded function. Then the function 
u^{x,y,t) = E/(Xf'^) satisfies the backward Kolmogorov equation associated with the 
SDE [?, Ch. 6] 

du'^ 1 

= {-Co + Ci)u' , with u'^\t=o = /• 



dt 



(A.2) 



We look for a solution of ( A.2| in the form of a power series in r: 

u^{x,y,t) = uo{x,y,t) + TUi{x,y,t) +T^U2{x,y,t) + 



(A.3) 



We substitute (A.3 1 in (A.2 1 and by equating equal powers in r we obtain the following 



sequence of equations: 
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= 0, 






C0U2 


= Ciui 



duo 
dt ' 
dui 

~dt' 



(A.4a) 
(A.4b) 

(A.4c) 



From (A. 4 1), since the process generated by £0 is ergodic we deduce that the first term in 
the expansion is independent of y, so uo{x,y,t) = uo{x,t) . In order for equation (A.4d) 
to have a solution we need 



and since uq is independent of y 



so that 



dug 
dt 

dt 



{y ■ ^xUo)p 

V{x) ■ VjjUq. 



This imphes that, to leading order, the dynamics are deterministic since this is the 
Liouville equation corresponding to the ODE x — v(x) as expected. We now calculate 
the first order correction in the expansion (A. 3). Equation (A.4d) becomes: 

- CqUi = y -VxUq - v{x) -V^cUq. (A. 5) 

We can solve this by setting 

""1 = x{x, y) ■ ^xUo + '^{x, t), 
where the ^! term belongs to the null space of £0 • We impose the normalization 

{X{x,y))p^0. (A.6) 
^.5 1 we obtain the cell problem 

(A.7) 



If we now substitute the expression for ui in (A. 5 1 we obtain the cell problem 

-^ox = y-v{x) 



which under the condition (A.6 1 gives the solution x — V ~ v{x) and so 

Ml = (2/ - ■ VxUq + ^{x, t). 

If we now substitute this expression for ui into ( A.4| :) we obtain 

-C0U2 y{y - v{x))'^ : V^V^uq - (yxv{x)y) ■ V^uq + y ■ 

and by applying the solvability condition we rnd up with 
' dui\ , , , y 



dui 
~dt 



dt 



{y{y ~ v{x)) : Vx^xUo - {\7xv{x)y) ■ VxUq + y ■ Vx^)p- 



Thus 



5* 



({y - v{x)){y - : V^V^uq - (yxv{x)y) ■ V^Uq + y ■ Va:^)p 

{v{x){y - v{x))'^ : W^W^Uo)p, 

-\/xv{x)v{x) ■ \/xUo + v{x) ■ Va;* + ^A^Mq. 
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TUi) we find that u = uq + t^! and that u satisfies 



v{x) ■ V^u-t{V^v{x))v{x) ■ V^u+ -A^u + 0{t^). 



(A., 



If we now set u — {uq 
du 

Ik 

This implies that the effective equation describing the corrections to the Langrangian 
dynamics is 

X = v{x) — t{'\/ xv{x))v{x) + y/rW , 
since (A. 8 1 is the backward Kolmogorov equation for this SDE. We reiterate that the 
derivation of equation (A. 8 1 is just formal and is not a proof that solutions to (4.10) are 
indeed close to solutions to ( A.l I. To prove such a result would require more sophisticated 
techniques such as those in Part III of [23) . 



Appendix B In this Appendix we present the two Lemmas needed for the proof of 

Theorem m _ 

Lemma 12 Let (j)(x,i/, At), (l){x,At) defined in equations (2.111, (2.201. Then there exist 
constants M, K, Ci independent of r, At, such that 

Py(j){x, y. At) < M and 
I \PJix, y, At) - cp{x, Ai) 1 1 < (1 + KAt) \\x-x\\ + Ci^ 

Proof: We start by analyzing the dterministic subequations ( |2.12[ ), in particular the 
equation for y. More specifically, we show that Py(j)j{x,y, A) =0(1), uniformly in r. We 
have 



Py(j)j{x,y,At) = yexp 
where 



(n+ l)r 



At 



exp 



At- s 
{n+l)T 



f,is)ds, (B.l) 



fj{t) = djVj a — b{n + l)r exp 



(n+ 1)t 



If we make the substitution q — e<"+i)^ then equation (B.l I becomes 



Py(j)j{x,y,At) = ye + (n + l)\/fe 

and now since v-j is bounded we deduce that 



dj Vj \ a 



b{n + l)r 



dq 



Py4).j{x,y,At) <ye <"+i)- + C(?i + 1) I 1 - e <"+i)- < Ciy/r 



We now study the a;-equation which can be written, using (2.201, as 
P^^j{x,y, At) ^x + V^y{l - eV^) + 







r At-s "1 
1 - e 


Jo 





fj{s)ds 



Note now that the integral on the right hand side of the above equation can be written 
as 



At 



1 — exp — 



At-s 
(n+ 1)t 



f,{s)ds 



At 



fj{s)ds 



At 



exp 



At-s 
(n+ 1)t 



fjis)ds 
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Now notice that the second part of the integral is the same that appeared in the y- 
equation multiphed by -y/r. Thus we obtain the bound 



At 



< Cir. 



(B.2) 



Now we have that 
^At 







fj{{ej,x) + (n+ l)r(ej,y)(l - e -))rfs = Atfji{ej,x)) 

/■At 

+/'((e„x)) / Ci(n + l)r(e„y)(l-e-^)ds 
Jo 

where we have taken a Taylor expansion around (ej,a;J(0)). Thus we have 
P,0, (S, y, At) = q},{x, At) + ry(l - e(^) 

/.At 

+ /'((e,,J)) / Ci{n+l)T{ej,y){l-e-^)ds 



(B.3) 



Using now the fact that Vj is Lipschitz continuous with constant Lj < L for all j and 
also bounded is easy to see that 

\\P^^j{x,y,At) - 4>j{x,At)\\ < (1 + LjAt)\\x-x\\ + CiV^ + C2TAt + Cgr^At (B.4) 



The proof of the first statement of the lemma now follows from ( B.2 1 , the bound \Py4>j At) \ < 
C, together with the fact that and Py(f){x,y, At) is given as a composition of every 



Py(t)j{x, y, At). In order to prove the second statement we need to use ( |B.4[ ) 

I \Pj{x, y, At) - Ai) 1 1 = I |P,0„ o ... 01 (J, y, At) - o . . . 01 (x, At) \ \ 

<(l + LAt)"P-a;|| + Ci^, 
<{l + KAt)\\x~x\\ + CxV^. □ 

Now we study the effect of the additive noise. We have the following lemma: 



Lemma 13 Le< 7, (7(7, ^, denote the random variables that we add in equations (2.11 1,(2.20 I. 
Then there exist a constant M independent of At, r 

nWPMl.Lt) ~ <jW{t)\\)^ < Mr (B.5) 



Proof: We can solve (2.161 to obtain 



c^(t) = 2;^(0) + ^/7(l-e(^^)?/(0) + c^ / (l - e T^)dW, 



Thus 



2;(t) = 2;(0)e('.+i)- + 



P. 



■9{l,i.t)^a (l-e \dW, 
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where 

Thus, and upon noticing that \/(i)7 = W{t), 



n\\p^9{i,i,t)-<jw{t)\\f=<j^ [ 

Jo 



lo 

2/ 



— — ^(1 - e ) < Mr, 



and the proof is complete. □ 
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